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ABSTRACT 

Context. The first observational evidence of a connection between supernovae (SNe) and gamma-ray bursts (GRBs) was found about 
a decade ago. Since then, only half a dozen spectroscopically confirmed associations have been discovered and XRF 1003 16D asso- 
ciated with the type-Ic SN 2010bh is among the latest. 

Aims. We constrain the progenitor radius, the host-galaxy extinction, and the physical parameters of the explosion of XRF 
100316D/SN 2010bh at z = 0.059. We study the SN brightness and colours in the context of GRB-SNe. 

Methods. We began observations with the Gamma-Ray burst Optical and Near-infrared Detector (GROND) 12 hours after the GRB 
trigger and continued until 80 days after the burst. GROND provided excellent photometric data in six filter bands covering a wave- 
length range from approximately 350 to 1800 nm. significantly expanding the pre-existing data set for this event. Combining GROND 
and Swift data, the early broad-band spectral energy distribution (SED) is modelled with a blackbody and afterglow component atten- 
uated by dust and gas absorption. The temperature and radius evolution of the thermal component are analysed and combined with 
earlier measurements available from the literature. Templates of SN 1998bw are fitted to the SN itself to directly compare the light- 
curve properties. Finally, a two-component parametrised model is fitted to the quasi-bolometric light curve, which delivers physical 
parameters of the explosion. 

Results. The best-fit models to the broad-band SEDs imply moderate reddening along the line of sight through the host galaxy 
( A t '.host = 1.2 ± 0.1 mag). Furthermore, the parameters of the blackbody component reveal a cooling envelope at an apparent initial 
radius of 7 x 10” cm, which is compatible with a dense wind surrounding a Wolf-Rayet star. A multicolour comparison shows that SN 
2010bh is 60-70% as bright as SN 1998bw. Reaching maximum brightness at 8-9 days after the burst in the blue bands. SN 2010bh 
proves to be the most rapidly evolving GRB-SNe to date. Modelling of the quasi-bolometric light curve yields M Ni = 0.21 ± 0.03 M Q 
and M C j = 2.6 ± 0.2 M a , typical of values within the GRB-SN population. The kinetic energy is £ t = (2.4 ± 0.7) x 10 52 erg, which is 
making this SN the second most energetic GRB-SN after SN 1998bw. 

Conclusions. This supernova has one of the earliest peaks ever recorded and thereafter fades more rapidly than other GRB-SNe, 
hypemovae, or typical type-Ic SNe. This could be explained by a thin envelope expanding at very high velocities, which is therefore 
unable to retain the y-rays that would prolong the duration of the SN event. 

Key words, gamma-ray burst:individual:XRF 100316D - supemovae:individuaI:SN 2010bh 


1. Introduction 

The core collapse of massive stars is thought to give rise to su- 
pemovae (SNe) and long-duration gamma-ray bursts (GRBs, du- 
rations larger than 2 s: Kouveliotou et al. 1993). The first clue 
to the connection between these events is the similarity in their 
kinetic-energy scale (see Woosley & Bloom 2006 for a review). 
Even before the discoveries of the late 90’s, a few authors dis- 
cussed the possible association between high-energy outbursts 
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and SNe (Colgate 1968; Paczynski 1986). In theory, the collaps- 
ing core of a very massive star, whose envelope has been blown 
away by its own stellar winds, i.e.. a Wolf-Rayet (WR) star, 
can lead to the formation of a relativistic jet that will produce 
high-energy emission (Woosley 1993; Woosley & MacFadyen 
1999) in the form of a GRB or an X-ray flash (XRF, softer 
events thought to be mainly produced by the same mechanism 
as GRBs; Heise et al, 200 1 ; Kippen et al. 2004; Sakamoto et al. 
2005; Heise & in 't Zand 2001). When the jet collides with the 
circumstellar material, it produces a multi-wavelength afterglow 
(for a review, see Zhang & Meszaros 2004). Theoretically, the 
energy of the core collapse should also be capable of expand- 
ing the remaining envelope. However, it is unclear how this can 
happen, or even if there is always enough energy for a SN ex- 
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plosion. The kind of SN thought to be associated with GRBs are 
those called stripped-envelope (SE) SNe, whose hydrogen enve- 
lope have mostly been removed. Events such as these are unusu- 
ally energetic stellar explosions and sometimes labelled “hyper- 
novae” in the literature (e.g., Paczynski 1998: Hansen 1999). 

To date, only type-Ic SNe have been related to GRBs. The 
first and most representative case was that of SN 1998bw, 
which was found to be associated with the soft GRB 980425 
(e.g., Galamaetal. 1998b; Kippenetal. 1998). A couple of 
days after the trigger. SN 1998bw was discovered (Galama et al. 
1998a; Sadler etal. 1998) inside the 8' error circle of GRB 
980425 (Soffitta et al. 1998) in the underluminous late-type 
galaxy ES0184-G82 (z = 0.0085; Tinney et al. 1998). Although 
initially controversial (Galamaetal. 1998b; Pianetal. 1998), 
the physical association between these objects was supported on 
temporal and spatial grounds by the slowly variable X-ray source 
at the position of the SN (Plan et al. 2000; Kouveliotou et al. 
2004). As a result of this connection, both the GRB and SN re- 
search fields, which until then had both evolved more or less 
independently, were revolutionised by a single event. 

Five years after SN 1998bw, the association between GRB 
030329 and SN 2003dh came to light through a clear spectro- 
scopic identification (Hjorth et al. 2003; Kawabata et al. 2003: 
Stanek et al. 2003; Matheson et al. 2003) and became the first 
truly solid piece of evidence in favour of the GRB-SN connec- 
tion. Moreover, after 1998 there have been many other spec- 
troscopic associations, such as the cases of GRBs 02121 1 (SN 
20021t; Della Valle et al. 2003), 020903 (Soderberg et al. 2005: 
Bersier et al. 2006), 031203 (SN 20031w; Malesani et al. 2004), 
050525A (SN 2005nc; Della Valle et al. 2006b), 060218 (SN 
2006aj; Pian et al. 2006), 081 007 (SN 2008hw; Della Valle et al. 
2008), 09 1 1 27 (SN 2009nz; Berger et al . 20 1 1 ; Cobb et al . 20 1 0; 
Filgasetal. 2011), and 101219B (SN 2010ma; Sparreetal. 
201 1) based on hypernova features in their spectra. 

Furthermore, late-time bumps in the light curves of GRB 
afterglows have been interpreted as SN signals, e.g., GRBs 
970228 (Reichart 1999: Galamaetal. 2000; Reichartet al. 
2000), 980326 (Castro-Tirado & Gorosabel 1999: Bloom et al. 
1999), 011121 (Bloom et al. 2002; Greiner et al. 2003), 020405 
(Price et al. 2003; Masetti et al. 2003), 040924 (Soderberg et al. 
2006; Wiersema et al. 2008), 041006 (Stanek etal. 2005; 
Soderberg et al. 2006). 050824 (Sollerman et al. 2007). 060729, 
and 090618 (both in Cano et al. 201 la) to mention a few. These 
bumps show consistency in terms of colour, timing, and bright- 
ness with those expected for the GRB-SN population, but they 
are usually faint, which hampers the spectroscopic identification. 
These re-brightenings have been detected in GRB light curves 
out to redshifts of ~ 1 (Masetti et al. 2005; Della Valle et al. 
2003: Bloom et al. 2009; Tanvir et al. 2010) owing to the sensi- 
tivity of current ground-based telescopes dedicated to follow-up 
observations. Sample studies of GRB-SNe (including bumps not 
spectroscopically identified) have been carried out to determine 
the luminosity distribution, the morphology of the light curves, 
and the physical parameters of the explosion such as kinetic en- 
ergy, ejected mass, and M \ i mass (Zeh et al. 2004; Ferrero et al. 
2006: Richardson 2009; Thone et al. 2011). It has been asserted 
that GRB-SNe are in general brighter than the local sample of SE 
SNe. except for cases such as SN 2010ay (Sanders et al. 2011). 


In contrast to the claims of Stanek et al. (2005), no correlation 
has been found between the brightness at maximum with the 
shape of the light curve (Ferrero et al. 2006). 

Additional information about the explosion can be obtained 
from its early emission. One particular case is that of the 
soft XRF 060218 associated with SN 2006aj (Campana et al. 
2006; Pianetal. 2006; Ferrero etal. 2006; Cobb etal. 2006; 
Modjaz et al. 2006; Sollerman et al. 2006; Mirabal et al. 2006). 
This displayed an early X-ray and UV emission, which 
was interpreted as thermal radiation produced by the shock 
breakout from the surface of the progenitor (Colgate 1974; 
Falk 1978; Klein & Chevalier 1978; Matzner & McKee 1999; 
Waxmanetal. 2007; Nakar&Sari 2010). The envelope is 
heated up and owing to expansion it starts to adiabatically cool, 
which then shifts the emission to UV and optical wavelengths. 
From the analysis of this a signal, it is possible to constrain 
the afterglow component, derive both the temperature and lu- 
minosity of the thermal component, and compute the appar- 
ent radius of emission (e.g., Thone etal. 2011). Other exam- 
ples include SNe without detectable -y-ray emission that never- 
theless exhibit adiabatic cooling in the UV/optical and/or X-ray 
observations: SN 2008D (Soderberg et al. 2008; Malesani et al. 
2009; Modjaz etal. 2009; Mazzalietal. 2008), SN 2008ax 
(Romingetal. 2009), SNLS-04D2dc (Schawinski et al. 2008), 
and SN 2010aq (Gezari et al. 2010). to mention a few. As addi- 
tional information, it has been claimed that high-energy emis- 
sion in GRB-SNe comes from accelerated shock-breakout pho- 
tons rather than highly relativistic jets (Wang et al. 2007). 

Whilst no SN signature is expected for short GRBs, which 
are thought to be produced by the mergers of compact objects 
(for deep non-detections, see the cases of GRB 050509B in 
Hjorth et al. 2005a and Bloom et al. 2006, and GRB 050709 in 
Fox et al. 2005 and Hjorth et al. 2005b), there are a few sup- 
posedly long events where an expected SN appearance was 
never detected. In the cases of GRBs 060505 (Fynbo et al. 2006; 
Ofek et al. 2007) and 060614 (Fynbo et al. 2006; Gal-Yam et al. 
2006; Della Valle et al. 2006a), there are very tight constraints 
on the SN signature, which go down to 1% as bright as SN 
1998bw. The validity of the non-detections of SNe components 
for these two GRBs as a classification tool is a point of con- 
troversy (e.g., Zhang et al. 2009). We refer to Kann et al. (2011) 
for an exhaustive discussion on GRBs 060505 and 060614, in 
the light of the optical luminosity of their afterglows. 

We analysed the optical and near-infrared (NIR) data of XRF 
100316D and its associated SN 2010bh. The paper is organised 
by summarising the observations, data acquisition, reduction, 
and analysis in Sect. 2. The main results are treated separately as 
three different sections. The modelling of the early broad-band 
spectral energy distribution (SED) provides the progenitor radius 
and the host-galaxy extinction and is presented in Sect. 3. The 
multi- wavelength light and colour curves are analysed in Sect. 4 
along with comparisons with previous GRB-SN events. In Sect. 
5, the quasi-bolometric light curve is analysed and the physical 
parameters of the explosion are derived. Finally we gather our 
conclusions in Sect. 6. 
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Fig. 1. HST/WFC3 FA\AW, GROND /'-band, and GROND host-subtracted /'-band images of the GRB-SN field, respectively, from left to right. 
Images were taken at around maximum brightness. Each panel is approximately 13" x 13" in size. The HST image shows significant galaxy 
structure near the explosion site (marked with two lines), which is blended with the object of interest in the GROND image shown in the middle 
panel. GROND images in the middle and right panels are shown using the same flux scale. 


2. Observations 

Gamma-ray emission from XRF 100316D triggered the Burst 
Alert Telescope (BAT; Barthelmy et al. 2005) on board the 
Swift satellite (Gehrels et al. 2004) on March 16, 2010, at to = 
12:44:50 UT (Stamatikos et al. 2010). It turned out to show a 
soft y-ray spectrum (Sakamoto et al. 2010) and a duration of at 
least 1300 s, one of the longest ever measured (Fan et al. 2011; 
Starling et al. 2011). About 1 5 h thereafter, a spectroscopic red- 
shift of 0.059 was published for the host galaxy (Vergani et al. 
2010a,b). Observations of the Gamma-Ray burst Optical and 
Near-infrared Detector (GROND; Greiner et al. 2007, 2008) 
confirmed that the new source became evident about 16 h af- 
ter the burst (Afonso etal. 2010). The rising of the supernova 
was verified photometrically by Wiersema et al. (2010) only 
three days after the trigger. The spectroscopic confirmation of 
Chornock et al. (2010b) came approximately six days after the 
burst, which was confirmed two days later by Bufano et al. 
(2010a). The SN was officially named SN 2010dh eight days af- 
ter the trigger (Bufano et al. 2010b: Chomock et al. 2010c). On 
March 26. 2010, additional GROND observations of SN 2010bh 
were reported along with results at the first attempts of host- 
galaxy subtraction (Rau et al. 2010). 

2.1. GROND 

The multi-channel imager GROND (Greiner et al. 2007, 2008), 
mounted at the MPG/ESO 2.2m telescope on La Silla, Chile, 
started observations of XRF 100316D 1 1.7 h after the trigger, 
simultaneously in g’r'i'z'JHK s . with an average seeing of 1 1 . 
as soon as the astronomical night began. 

The GROND data were reduced using standard 
pyraf/IRAF 1 tasks (Tody 1993), similar to the procedure 
described in detail by Kriihleretal. (2008). A general point- 

1 IRAF. the Image Reduction and Analysis Facility, is distributed by 

the National Optical Astronomy Observatory, which is operated by the 
Association of Universities for Research in Astronomy (AURA), Inc., 
under cooperative agreement with the National Science Foundation 
(NSF); see http://iraf.noao.edu. 


spread function (PSF) model was constructed using bright field 
stars, from which the full width at half maximum (FWHM) was 
derived for each image. Aperture photometry was performed 
for science and calibration objects using an aperture size equal 
to the FWHM. The SDSS field at coordinates R.A.(J2000) 
= 06 h 59 m 33 s .6, Dec.(J2000) = -17°27'00" was observed 
during photometric conditions to calibrate our images by 
performing relative photometry. A total of six stars in the field 
of XRF 100316D were employed for this purpose (see Table 
A.l of the online material). The same set of stars was used 
to calibrate JHK against the 2MASS catalogue. Calibration 
uncertainties vary in the range 0.002 - 0.020 mag for g'r'i'z' 
and 0.02-0.12 mag for JHK 5 , which together with catalogue 
systematics are added in quadrature to the statistical error. 

A deep host-galaxy observation was carried out on 
November 5, 2010. This observation resulted in images with 
mean seeing of 0774 and was used as a reference image for 
the subtraction of the host contribution from the early epochs. 
To develop a notion of the host contamination, we refer to 
Fig. 1 , where the position of the transient relative to the host 
galaxy is shown along with an F& 14 W image of the Wield 
Field Camera 3 (WFC3) on board the Hubble Space Telescope 
(HST) as reference. In addition, the coordinates of the cen- 
tre of the host galaxy, catalogue name Anon J07 103 1-56 15, 
were obtained from the last GROND observation: R.A.(J2000) 
= 7 h 10 m 30 s .37(±0 s .07), Dec.(J2000) = -56°15'2072(±073). 
To align the input and reference images, we use the WCSREMAP 
package 2 . For the main purpose of subtracting template from sci- 
ence images, the H0TPANTS package- 1 was employed. It matches 
the point spread function (PSF) and count flux of both input im- 
ages. It uses Gaussian functions to model the PSF in sub-regions 
of the original image. Since reduced images are essentially the 
combination of at least four g'r'i'z' and 24 JHK s dithered in- 
dividual exposures, three and five Gaussian functions of differ- 
ent widths and degrees of freedom are employed for g'r'i’z’ 
and JHK, respectively. Point sources define the transformation 

2 http://www.astro.washington.edu/users/becker/wcsremap.html 

3 http://www.astro.washington.edu/users/becker/hotpants.html 
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that is used to scale the template image to the flux of the sci- 
ence image. The routine outputs a noise map of the resulting 
difference image, which is employed to derive the uncertainties 
in the measured fluxes. The package DAOPHOT/IRAF was used 
to compute the magnitudes and their corresponding errors. This 
procedure was executed for a total of 140 individual images in a 
total of 20 epochs (see Fig. 1 for an example of image subtrac- 
tion in the i'band). The position of the transient is R.A.(J2000) 
= 7 h 10 m 30 s .55(±0\05), Dec.(J2000) = -56°15'20"0(±0"2) in 
host-subtracted optical images. The resulting photometry is tab- 
ulated in Table A. 2 of the online material. 

2.2. Swift /XRTand UVOT 

On board the Swift satellite, the X-Ray Telescope (XRT; 
Burrows et al. 2005) and the UVOT (Roming et al. 2005) started 
observations of XRF 1 003 1 6D at to + 2.4 min (Stamatikos et al. 
2010). Whilst a bright X-ray source was detected inside the BAT 
error circle, initially no afterglow candidate was found by UVOT 
(Oates et al. 2010). 

However, in deeper images taken at 63 + 33 ks in the uv w\ 
filter and at to + 63 ks in the u band, we found evidence of emis- 
sion in excess of the host-galaxy contribution. To remove the 
contribution from the host galaxy, we requested ToO observa- 
tions in the uvwl and u filters, which were taken at to + 3 x 10 7 s 
(347 d after the burst), and amounted to a total exposure time of 
1525 and 1369 s, respectively. We measured a host galaxy con- 
tribution within the source aperture of 25 ± 2 and 39 ± 3 giy 
in the uvwl and u bands, respectively. Subtracting this contribu- 
tion from our earlier-time data gave us a 3cr detection in 1 894 s 
of uvwl- and in 4901 s of u-band data taken at mid-times of 
to + 33 ks and to + 63 ks respectively. In our analysis, we in- 
clude only the uvwl detection at to + 33 ks, since there were no 
GROND data contemporaneous with the epoch of our u-band 
detection. All UVOT data were reduced following the procedure 
described in Poole et al. (2008). To minimise the contamination 
from the underlying host galaxy, the source flux was measured 
within a circular source-extraction region of a 375 radius. An 
aperture correction was then applied in order to remain compat- 
ible with the UVOT effective area calibrations, which are based 
on 5" aperture photometry (Poole et al. 2008). The background 
was taken from a source-free region with a 15" radius close to 
our source. 

The relatively bright X-ray afterglow (30-40 ent s ” 1 be- 
tween 144 and 737 s after the trigger) faded considerably at the 
beginning of the second XRT epoch(4)+33 ks; see Starling et al. 
20 1 1 . for a detailed analysis). In the subsequent analysis, we em- 
ployed XRT data at stages contemporaneous to GROND obser- 
vations. specifically in the interval from 33 to 508 ks after the 
burst. These data were obtained from the public Swift archive 
and reduced in the standard manner using the xrtpipeline 
task from the HEAsoft package, with response matrices from 
the most recent CALDB release. Ali data were obtained in photon 
counting mode and downloaded from the XRT light curve repos- 
itory (Evans et al. 2007. 2009). Spectra were grouped using the 
grppha task. 

All the data discussed throughout the paper were corrected 
for the Galactic-foreground extinction. of E(B - F)oai = 0.117 


mag with Ry = 3.08 (Schlegel et al. 1998). All uncertainties in 
the following analysis are quoted at the ler confidence level. 

3. Early broad-band SED 

Using data between 33 and 54 ks (roughly from 11 to 15 h) af- 
ter the burst, we compiled two early broad-band SEDs of XRF 
100316D. GROND provides detections in g'r'i'z'J and up- 
per limits in HK S , whilst the count rate of contemporaneous 
Swift/XKT observations is already 0.01 ent s" 1 at to + 33 ks and 
decaying. After combining data from 33 to 508 ks after the trig- 
ger. XRT provides only three bins in the X-ray energy range. As 
follows, XRT data are scaled to the GROND first two epochs by 
using a decay index of a = -1.3 ± 0.2, which is derived from 
the same XRT data over time. 5tw/(/UVOT observed only in the 
uvwl filter at these stages. 

3.1. Modelling scheme 

Early blue emission coming from the XRF position was detected 
by GROND. For data acquired 42.5 ks after the burst, we mea- 
sured a colour g' - r' = -0.30 ± 0.06 mag, in contrast to the 
red afterglows that usually follow a GRB (photon index F in the 
range 1. 2-2.5). Similar observations were made by Cano et al. 
(201 lb), who found that this emission is incompatible with syn- 
chrotron radiation. The adiabatic cooling of the expanding atmo- 
sphere following the shock breakout gives us a reasonable expla- 
nation of the observed blue colours (e.g., Cano et al. 201 lb). In 
this scenario, the emission from the shock breakout lasts only 
a few hours after the core collapse and its SED resembles a 
blackbody at a high temperature of the order of 10 6 K (or 0.1 
keV equivalently; e.g., Campanaetal. 2006). Thus, we inter- 
pret the blue colours in our observations as the thermal compo- 
nent associated with the cooling of the shock breakout, which is 
similar to what has been claimed for XRF 060218/SN 2006aj 
(Campana et al. 2006; Waxman et al. 2007) and other early- 
caught SNe (e.g., Soderberg et al. 2008; Modjaz et al. 2009 on 
SN 2008D; Roming et al. 2009 on SN 2008ax). In contrast to the 
idea of a thermal component producing the observed emission, 
there are no significant contemporaneous detections of UVOT 
in the uvwl bandpass, which alludes to high reddening. On the 
other hand, i'-J= 0.06 ± 0. 1 5 mag measured at the same epoch 
appears to have an underlying additive red component, which is 
represented in this case by the GRB afterglow synchrotron emis- 
sion (see Fig. 2). 

To test this hypothesis, we modelled the SEDs at 42.5 and 
50.0 ks (the latter shown in Fig. 2) using two additive compo- 
nents: (1) a power law for the afterglow, and (2) an ideal black- 
body for the thermal component, respectively, of the forms 

PeTT) = C\ E~ l and (1) 

CzE-dE 

B E (T m ) - {kTm )4{e£n+zU/BF _ j 

where C\ and Cj are normalisations, E is the spectral energy 
in units of keV, F is the photon index, k is the Boltzmann con- 
stant, 7 bb the intrinsic blackbody temperature, and z = 0.059 is 
the redshift (Chornock et al. 2010c). The model also accounts 
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Table 1. Fits to broad-band spectra constructed using GROND and Swift/XKT data. 


Time Interval 

Reddening 

E(B- T) host 

L7bb 

Lbb 

Rbb 

r 

'Vh+osi 

x 2 /r 

[ s] after the trigger 

Law 

[magj 

[eV] 

[ 10 47 ergs'” 1 ] 

f 10 12 cm] 


[ 10 22 cm- 2 ] 


42182-42879 

MW 

0 21 + ? J1 
-0.21 

81!$ 

4-l!|? 

27±8 

1.82+0.05 

4.2+0.5 

4.9/4 

42182-42879 

LMC 

0 19 +0I ° 
u ‘ -0.19 

81!$ 

4.°! 1 : $ 

26±8 

1.82+0.05 

4.2+0.5 

4.9/4 

42182-42879 

SMC 

0 16 +on 

-0.16 

82!$ 

3.4 ± 1.7 

T 4+ 7 
^ -6 

1.82+0.05 

4. 1+0.5 

4.9/4 

42182-42879 

MW 

0.39 fixed 

79 ±3 

7.4 ± 0.9 

39 ±4 

1.73 + 0.06 

4-4!“ 

8.7/5 

46630-53807 

MW 

0.39 ± 0.03 

78 ±2 

7.2 ± 1.1 

38 ±3 

1.77 + 0.05 

4.4 ± 0.4 

13/5 

46630-53807 

LMC 

0.38 ± 0.03 

78 ± 2 

7.0 ± 1.0 

38 ±3 

1.76 + 0.06 

4.4 ± 0.4 

16/5 

46630-53807 

SMC 

0.37 ± 0.04 

79 ±2 

6.2 ± 0.9 

35 ±3 

1.77 + 0.05 

4.4 ± 0.4 

18/5 


Notes. The model consists of a blackbody plus a power law, both attenuated by optical/NIR reddening and X-ray metal absorption. 
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Fig. 2. Broad-band SED at 50 ks after trigger. The observed data are 
represented by black filled circles. The HK S bands provide only 3cr 
upper limits shown as arrows. The thick grey line shows the best-fit 
model able to reproduce the data: extinguished power law plus black- 
body components, which are shown individually with dotted and dashed 
grey lines, respectively. The thin continuous grey lines show the unex- 
tinguished versions of the power law corresponding to the afterglow and 
the thermal component. In the lower panel, the residuals of the best fit 
are plotted. 


for host-galaxy extinction based on either Milky Way (MW; 
Ry = 3.08). Large Magellanic Cloud (LMC; Ry = 3.16). or 
Small Magellanic Cloud (SMC; Ry = 2.98) extinction laws and 
soft X-ray metal absorption. The Galactic metal absorption is 
fixed to be A^n.Gai = 7.05 x 10 20 cm -2 (Kalberla et al. 2005). 
The model has a total of six free parameters: C i, Ci, F. T m . 
E(B - T)hosi, and /Vr host- The luminosity of the blackbody is 
computed as Lbb = 8.0525 Cz(l + z)D~ 0 erg s~ ! , where D V: is 
the luminosity distance to the transient in units of 10 kpc. A lu- 
minosity distance of 240 ± 17 Mpc to the host galaxy of XRF 


100316D is employed as computed by the NED database 4 fol- 
lowing the standard ACDM model, using the redshift measured 
by Chomocketal. (2010c), a Hubble constant of 74.2 km s ' 1 
Mpc ” 1 (Riessetal. 2009), and the model of the local velocity 
field by Mould et al. (2000). 

The different set of parameters are summarised in Table 1. 
All best-fit parameters are consistent between the two epochs 
within their statistical uncertainties between both epochs. The 
temperature 7 bb and luminosity Lbb of the thermal component 
are in the range 78-81 eV and 4 - 7 x 10 47 erg s _l from Cols. 4 
and 5 of Table 1 . respectively. The afterglow power-law pho- 
ton index F 1.8 is shown in Col. 7. The reduced ^ 2 (or 
X, A 2 = x 2 Ie- where (t is the number of degrees of freedom) im- 
proves from 16.6 for the model without extinction to 2.7 for the 
model extinguished by MW-like dust. The best fit to the sec- 
ond SED epoch is shown in Fig. 2 with a thick grey line. The 
UV dust feature characteristic of the MW extinction law gives 
the more precise results, although it is poorly constrained at the 
bluer end. The X-ray tail of the blackbody fits the two data bins 
at around 1 keV, whilst the power law fits the only data point 
at * 3 keV; both components contain significant absorption. We 
note that the fitting of the X-ray data has practically no residuals, 
i.e., the model over-predicts the data in this energy range. 

In addition to the luminosity and temperature, it is possible to 
compute an apparent emission radius of the thermal component 
(L?bb) from these two measured quantities by assuming isotropic 
radiation from an ideal blackbody. For all our trial models, the 
corresponding radii were calculated as /?bb = ( 47 r<rL^ B r B 4 B )~ 1/2 , 
where <x is the Stefan- Boltzmann constant (Col. 6 of Table 1). 
Our best-fit parameters yield a radius of 3-4 x 10 13 cm. which 
is two orders of magnitude larger than the typical sizes of the 
most likely GRB progenitors (WR stars; Cappa et al. 2004). 

The amount of metal absorption (estimated based on an 
equivalent hydrogen column density at solar metallicity) and 
dust extinction required for a good fit are No. .host » 4 x 10 22 
cm " 2 and E(B- 1/ )h 0 st = 0.2 -0.4 mag, respectively. These two 
parameters are fitted independently and with no assumption be- 
ing made about the environment gas-to-dust ratio. However, the 
available data is not enough to allow us to distinguish among 
the different extinction laws employed in the modelling proce- 
dure. Whilst the reddening is consistent with values found in 

4 http://nedwww.ipac.caltech.edu/ 
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previous studies (Starling et al. 201 1; Cano et al. 201 lb), the hy- 
drogen column density differs significantly from the results of 
Starling et al. (201 1). The discrepancy is due to the additional 
constraints provided by the optical/NIR data, which were in- 5? 
eluded to derive the parameters of the blackbody contribution. A 
With our data set it is possible to tie the thermal component 
at both low and high energies, which provides a more accurate 
value for the absorption by heavy elements because of the greater 
constraint on the X-ray flux. 

3.2. Host-galaxy extinction 

We then attempted to evaluate the constraint on the extinction ” 
value derived in the previous section. Contour plots of redden- « 
ing against blackbody temperature, luminosity, and spectral in- *0 
dex from the best fit to the data are shown in Fig. 3, respectively. '~Z 
The ler and 2cr contours show that the solution to our best-fit J 
relation is well-defined in all three panels. In the first panel there 
is a second solution to our best-fit relation at the 3cr level, which 
is at lower temperatures and similar reddening, although its sig- 
nificance is rather low. In the middle panel, the contours show 
a slight trend of proportionality between luminosity and extinc- £ 
tion, which is expected. At luminosities higher than 1.1 x 10 48 J 
erg s~* , there is no possible solution, because the model becomes & 
much brighter than the X-ray data and (Vn.host cannot be allowed J 
to vary and compensate for the X-ray luminosity of the model. £ 
Contours of the variations in the reddening and spectral index are 
shown in the lower panel, where the extinction is again restricted 
to high values. 

Analyses of optical spectroscopy of the host galaxy have 
found that the line-flux ratio of Ho to H f3 in the Hu region 
coincident with the SN position is consistent with zero extinc- 
tion (Starling et al. 2011; Levesque et al. 2011). Nevertheless, 
the Aj/host ~ 0 estimate from the spatially resolved spectroscopy 
of Levesque et al. (201 1) probes a much larger region (« 1.3 
kpc 2 ) than the one probed by our line of sight to the GRB-SNe. 

A sufficiently high dust dumpiness could explain our high ex- 
tinction values along the line of sight and Ai/ hosl ~ 0 when in- 
tegrated over a larger patch. The position at which the hydrogen 
lines are formed, i.e., the H 11 region, might be located in front of 
the explosion site and not probe the same line of sight through 
the host galaxy. 

Alternatively a broken power-law model was fitted to the 
data in Fig. 2 for which Xu 2 = 3.0, which is slightly larger 
than that of the model that consists of a blackbody plus power 
law. The broken power law provided a closer fit when we as- 
sume that there is no host-galaxy extinction, significantly differ- 
ent from our results for the blackbody plus power-law model. 
However, the low-energy spectral slope of j3 = +0.5 is incom- 
patible with synchrotron radiation (/J ra ax = + 1 /3; e.g.. Sari et al. 

1998). Furthermore, the spectral break lies between the g'and 
the uvwl bands at Vbreak = (8 ± 3) x I0 5 GHz, which is incon- 
sistent with the self-absorption feature usually observed at radio 
frequencies (v a ~ 2- 13 GHz; Galama et al. 1998c; Taylor et al. 

1998; Granotetal. 1999; Galama etal. 2000). Given also that 
Atthost > 0-2 mag is derived in Sect. 3.1 and preferred by other 
authors (see next paragraph), the broken power- law model can 
be discarded with confidence. 



0 0.1 0.2 0.3 0.4 0.5 


E(B-V) 

Fig. 3. Contour plots for the best-fit parameters from the modelling of 
the second epoch. From the outer- to the innermost contours, the green, 
red, and black lines are 3, 2, and l<r contours, respectively. The tail at 
low F in the bottom panel is caused by numerical inaccuracies. 


Additional evidence of large reddening along the line of 
sight inside the host galaxy was found by Starling et al. (2011), 
£(B - F)host ~ 0-9 mag. The reddening was derived by fit- 
ting Swfft/BAT+XRT data and a u-band 3<x upper limit pro- 
vided by SwififUVOT data in the interval from 638 to 737 s after 
the burst. We employed a model that consists of a blackbody 
plus power law extinguished by SMC dust. A similar method 
was used here and by Campana et al. (2006), who determined 
E( B - F)host = 0.20 for SN 2006aj. Another attempt at dered- 
dening SN 2010bh was carried out by Cano et al. (201 lb), who 
found E(B- F)host = 0.18 ± 0.08 mag by assuming that the 
colours of SE SNe are all the same ten days after the F-band 
maximum brightness (Drout et al. 2010). Whilst this method is 
supported for the hydrogen atmospheres of type-IIP SNe by 
a line of physical arguments (e.g.. Olivares E. et al. 2010), it 
is entirely empirical for SE SNe. Nevertheless, their redden- 
ing value is larger than zero with a significance of 2.3cr and 
consistent with our calculations for our second-epoch SED of 
E( B~ F)host = 0.39 ± 0.03 mag at the 2.5<r confidence level. 

After correcting for E(B - F)host = 0.39 ± 0.03 mag. we 
obtained uvwl = 18.15 ± 0.27 and u = 19.70 ± 0.45 mag in 
the AB system at 33 and 63 ks after the burst, respectively. For 
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SN 2006aj (Campana et al. 2006). these values were uvwl = 
17.73 ± 0.21 and u = 17.77 ±0.15 mag at the same phase and 
redshift of SN 2010bh. If we had assumed that we should have 
seen a cooling envelope for SN 2010bh of comparable bright- 
ness and evolution as that shown by SN 2006aj, the host-galaxy 
extinction should have been higher than that estimated by our 
method. Nevertheless, the uncertainties are large and the com- 
parison of the uvwl measurements is consistent at the 1 ,2cr con- 
fidence level. Although in the u band the uncertainty is even 
larger, there is no consistency with the brightness of SN 2006aj 
at the 4cr confidence level and a higher host-galaxy extinction 
(±4 v'host ~ 2. 1 mag) would be required to reach the lcr level of 
consistency. 

In conclusion, we use the extinction value from the fit to 
the second-epoch broad-band SED throughout the paper, E(B- 
V )host = 0.39 ± 0.03 mag for MW-like dust with Ry = 3.08, 
given that compared to the first epoch the statistical errors in the 
g'r'i'z'J photometry are smaller and additional uvwl photom- 
etry is available. Moreover, when fixing the host-galaxy redden- 
ing to E(B - lOhos, = 0.39 mag, the quality of the fit to the 
first epoch is still acceptable (Xy 2 = 1-6), whilst the fit to the 
second epoch using E(B - I/) hos t = 0.21 mag from the first- 
epoch modelling results in^ 2 = 3.7 and unphysical parame- 
ters. Hence, values of host-galaxy extinction including the cor- 
rection for redshift (/(-correction based on the spectral model) 
for the GROND filters and their corresponding statistical un- 
certainty are A^ost = 1.50 ± 0.12, A r '.h 0S t = 1.10 ± 0.09. 
A'.host = 0.80 ± 0.08, A^host = 0.60 ± 0.06, A/j, os t = 0.39 ± 0.04. 
Ah , host = 0.22 ± 0.02, and A/r s ,h 0st = 0.14 ± 0.01, all in units of 
magnitude. 

3.3. Evolution of the thermal component 

Having determined temperature, luminosity, and radius for the 
thermal component, it was then possible to study their evolution. 
Starling et al. (2011) analysed combined BAT+XRT data until 
737 s after the trigger and derived blackbody temperatures, data 
that we used in the following analysis. Temperature and radius 
over time are shown in Fig. 4. The decrease in temperature and 
the increase in radius are both trends that are consistent with 
the cooling of the envelope due to expansion and, thus, with the 
explosive scenario. 

The temperature and radius of an expanding envelope that is 
cooling adiabatically have been theoretically shown to evolve as 
power laws (Waxman et al. 2007: Nakar & Sari 2010). We there- 
fore fit power laws to the decay and rise of temperature and ra- 
dius to the combined data presented here and in Starling et al. 
(201 1). The expression Tbb (O = T~Kt s was fitted to the evolv- 
ing temperature (solid line in the upper panel of Fig. 4). which 
gave a decay index of 5 = 0.3 + 0.2 for an initial temperature of 
Ti = 0. 1 7 ± 0.04 keV with 2 = 0.8. where k is the normalisa- 
tion of the power law. 

A model of the form /?bb(0 = Ro + vt y was employed to 
fit the radius measurements. In the case of linear expansion, i.e., 
y = 1. the radius grows at a mean velocity of v » 8.000 km 
s -1 between BAT+XRT (until 737 s) and GROND observations 
(at 42-54 ks after the burst). When y was allowed to vary, we 
obtained a growth index of y = 1.4 ± 0.3 and an initial radius 


CO 



Time after the GRB [s] 

Fig. 4. Temperature and radius evolution of the blackbody compo- 
nent. Black data points are those derived using combined BAT+XRT 
data (Starling et al. 2011) plus the temperature determinations using 
GROND data (« 46 ks). The first data bin from Starling et al. (2011) 
represents the interval from -175 to 144 s, however, due to fitting and 
plotting purposes, here it corresponds to 0- 144 s. Solid, dashed and 
dotted lines are different power-law models (see main text). Best fits 
to temperature and radius are shown in grey and brown solid lines, re- 
spectively. Grey crosses (upper panel) and open circles (lower panel)) 
are measurements of XRF 06021 8/SN 2006aj taken from Kaneko et al. 
(2007) and Campana et al. (2006), respectively. 


of Ro = (7.0 ± 0.9) x 10 11 cm with^ 2 = 1-1 (solid line in 
the lower panel of Fig. 4). The resulting radius is slightly larger 
than the size of WR stars (~ 10 1 1 cm; Cappa et al. 2004), which 
are thought to be the progenitors of long-duration gamma-ray 
bursts and type-Ic SNe (Woosley et al. 2002). Because of this, 
the initial emission radius of the thermal component might indi- 
cate the position at which a preexisting dense wind surrounding 
the progenitor becomes optically thin (e.g.. Campana et al. 2006: 
Soderberg et al. 2008; Balberg & Loeb 201 1). From the theoret- 
ical point of view, it is y = 0.8 (Waxman et al. 2007), although, 
the data do not favour this solution delivering^ 2 = 8 (dotted 
line in the lower panel of Fig. 4). 

For an assumption of adiabatic cooling, the luminosity must 
be constant and was fixed to Z.bb = 3.5 x 10 45 erg s -1 , value de- 
rived from early-time X-ray measurements ( Starling et al . 20 1 1 ) . 
Assuming the best-fit model for the radius evolution from the 
previous paragraph, we computed the evolution of the adiabatic 
temperature as Tbb = (47rcrZ^g R BB )~ 1/4 , which is shown by the 
dashed line in the upper panel of Fig. 4. This model is clearly 
inconsistent with our data set; however, it indeed shows consis- 
tency with the XRF 06021 8/SN 2006aj data set (grey crosses 
in the upper panel of Fig. 4; from Kaneko et al. 2007). Since 
the blackbody luminosity changes to 4-6 x 10 4 ' erg s _1 in our 
late-time measurements, the assumption of a constant luminos- 
ity is invalid and inconsistent with the late-time temperature de- 
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termination. This result implies that either the cooling is not 
strictly adiabatic for which energy injection from the inner core 
is needed, or 7 bb and Z-bb are overestimated by our modelling 
procedure, or underestimated by Starling et al. (201 1). 

4. Multicolour evolution of SN 2010bh 

We now present the entire GROND data set, which includes data 
in seven different bands covering the wavelength range from 380 
to 2300 nm. Photometry was corrected for host-galaxy extinc- 
tion computed at the end of Sect. 3.2. In addition, the after- 
glow component derived in Sect. 3.1 was subtracted from the 
data. Assuming a power-law decay for the optical afterglow with 
a = - 1 .3 ± 0.2 computed from the contemporaneous X-ray data, 
a flux and its corresponding uncertainty were derived from the 
power-law model fitted to the early SEDs for each epoch. We 
note however that typically ax £ «op!- The correction for the 
afterglow contribution is more significant at earlier times and in 
redder bands and it has a lower significance at around maximum 
brightness for all bands. 

4. 1. Optical and near-infrared light curves 

Figure 5 shows the light curves in the optical g'r'i'z' md the 
NIR JH bandpasses. The g'r’i'z' J light curves show the usual 
pattern of SE SNe: the redder the filter, the later and broader 
the peak. The H band peaks a few days earlier than J, although 
given the large uncertainties, the time difference is insignificant. 
The K s band shows no credible detections in any of our observa- 
tions down to limits in the range of 18.4- 18.9 mag (AB system) 
despite our observations covering the expected peak of the SN 
(10.5-25.5 d after the burst). In the following section, we anal- 
yse our data of SN 2010bh, one of the best-observed GRB-SNe 
to date, and compare the light curve with those of other SNe 
connected to GRBs. 

4.2. Colour evolution 

Six GROND colour curves are presented in Fig. 6 to analyse 
the colour evolution of SN 2010bh. Most colours show drastic 
evolution, where the cooling of the SN photosphere is evident as 
colours become redder. Non-significant variations over time are 
shown in r'- / 'and J-H colours. 

4.3. Comparison with SN 1998bw 

To study the luminosity evolution of SN 2010bh to other SNe. 
we fit SN 1998bw templates derived from the publicly avail- 
able UVBRI photometry (Galama et al. 1998b) to each of our 
GROND filters. The process of compiling these templates from 
the observed light curve of SN 1998bw is based on Zeh et al. 
(2004). Given their modus operandi, the NIR templates are inac- 
curate because they rely on an extrapolation of the t/KBB/data. 
Therefore, the three epochs of JHK s data in Patatet al. (2001) 
were used to define the zero points of the flux scale. For all fil- 
ters. it is assumed that the host-galaxy extinction of SN 1998bw 
is equal to zero (Patatet al. 2001; Clocchiatti et al. 2011). An 
analytical function was then employed to parametrise the SN 
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Fig. 5. Multi-colour light curves corrected for host-galaxy extinction 
and afterglow-subtracted. Filled circles represent detections and arrows 
are upper limits. Solid lines represent the overall fits and dashed lines 
individual components. Fits of SN 1998bw templates are extrapolated 
from a Sj + 60 d. For reasons of clarity, light curves were shifted along 
the magnitude axis and the systematical error in A y Jimt was not added 
to the error bars. 

1998bw templates. Two of the seven template parameters rep- 
resent the brightness and morphology of the light curve, the fac- 
tors of luminosity k and stretch s, respectively. These were de- 
fined provided that k = 1 and s = 1 for SN 1998bw (Zeh et al. 
2004). In addition, it was necessary to include a delay parameter 
ideiay that acts by shifting the template linearly in lime, i.e., as 
a time offset. The delay parameters that differ from zero do not 
necessarily mean that the onset of the SN is out of phase. 

Given that the early blue emission is inconsistent with the 
templates even after subtraction of the afterglow component de- 
rived in Sect. 3.1, an empirical power-law component was re- 
quired (dashed in Fig. 5). The power-law slope was fixed to 
a = -1.3. Leaving a free did not improve the fit and differ- 
ent values of a only negligibly affected the stretch and luminos- 
ity factors derived for the SN component. Table 2 summarises 
the results of fitting the SN 1998bw templates and the empirical 
power-law to the GROND data. The overall fits are shown in Fig. 
5 for each band using solid lines. 

Luminosity factors listed- in the first line of Table 2, which 
include the uncertainty in Ap.host, reflect differences between the 
colours of SN 2010bh and those of SN 1998bw. The g'band is 
as bright as SN 1998bw and, as a rough comparison, disagrees 
with the fainter B-band results from Cano et al. (201 lb) at the 
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Table 2. Fits of SN 1998bw templates to SN 2010bh. 



g' 

r' 

/" 

z f 

J 

H K s 

SN amplitude (k) a 
Stretch factor (s) 

Peak delay (/delay)* 

1.10 + 0.12 
0.77 + 0.01 
-3.3 

0.65 ± 0.05 
0.78 + 0.01 
-3.2 

0.54 ± 0.04 
0.65 ± 0.02 
-2.6 

1 .25 ± 0.07 
0.73+0.02 
-0.4 

0.69 ± 0.07 
0.84 ± 0.04 
-1.0 

0.63 + 0.08 <1.4 

0.67 + 0.07 
0 


3 For afterglow-subtracted and host-extinction-corrected data (/Svjioa = 1-2 ± 0.1 mag). 
b Relative to maximum luminosity of SN 1998bw after applying a time offset due to the stretch factor. 
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Fig. 6. Colour curves corrected for host-galaxy extinction and 
afterglow-subtracted. Filled circles represent detections in both filters 
and arrows are upper limits derived from detections in one of the two 
bands. Solid lines represent the colour evolution of SN 1998bw derived 
from the light-curve fits in Sect. 4.3. For clarity, error bars do not in- 
clude the systematic error in A^osi and for the last epoch, upper limits 
are not shown. 


5cr confidence level after including the uncertainty in their ex- 
tinction determination. Our r'i ' luminosity factors are also larger 
than previously reported from optical photometry of SN 2010bh 
(k ~ 0.4-0. 5: Cano et al. 201 lb) mainly owing to the use of 
a different extinction correction. There is certainly some intrin- 
sic blue excess in the g'band. although the difference from the 
r’ band of 40% has a rather large uncertainty of 13%. Higher 
luminosities at bluer wavelengths indicate that the photosphere 
has a higher temperature at its peak luminosity than that of SN 
1998bw. Given that the expanding atmospheres of SNe cool 
down over time and SN 20 lObh peaks about a week earlier than 
SN 1998bw in g\ it is a natural conclusion that it must be hotter 


at maximum brightness. It is also possible that a g'-band peak 
as early as 8 d after the trigger still contains a non-negligible 
contribution from the cooling shock breakout. We also note that 
the z '-band luminosity factor is much larger than in the other 
bands; it is brighter than SN 1998bw at the 3<r confidence level. 
After inspection of the spectra around maximum light presented 
by Bufano et al. (201 1), the Can 48579 emission line was found 
to contribute roughly 2-4% to the continuum flux integrated in 
the z'-band sensitivity range (8254-9528 A). The contribution 
of this spectral feature is not enough to account for the excess 
of 45% in the z'band relative to the infrared. The discrepancy is 
therefore attributed to the systematical uncertainties introduced 
when extrapolating the UBVRI data from Galama et al. (1998b) 
to construct the templates of SN 1 998bw. In summary, the lumi- 
nosity of SN 2010bh is a factor 0.5 -0.7 fainter than SN 1998bw 
in optical r'i' bands and 0. 6-0.7 in NIR JH bands. We also 
note that at late times SN 2010bh fades more rapidly in the op- 
tical than SN 1998bw did. although the fluxes if the templates 
were extrapolated after ~ to + 60 d. 

The stretch factors listed in the second line of Table 2 range 
from 0.6 to 0.8 and are at the low end of the GRB-SNe distribu- 
tion (see figure 5 of Ferrero et al. 2006). These are also consis- 
tent with the findings by Cano et al. (201 lb) for SN 2010bh. The 
optical light curves are wider than the stretch factor suggests, or 
in other words, the stretch factor predicts later peak times for the 
optical light curves. To account for this, tj e lay shifts the templates 
to earlier times by about three days in g'r'i 1 , which means that 
the optical light curves of SN 2010bh peak even earlier than the 
stretch factor suggests. After the inclusion of the delay parame- 
ter /del ay * the stretch factors now solely represent the width of the 
light curves and not the time of the maximum brightness com- 
pared to SN 1998bw. Moreover, peak times of roughly 8-9 d in 
the g'r ' bands (equivalent to the V in the Johnson filter system) 
corresponds to the earliest and fastest light curves of GRB-SNe 
observed to date (see table 2 of Richardson 2009). 

Offsets in peak time, in our case represented by feiay, can 
be explained as the result of delayed black-hole formation 
(Vietri & Stella 1999). In this scenario, the XRF might be trig- 
gered by the core collapse of the progenitor to a neutron star, 
soon after which accretion holds. The supernova would then oc- 
cur after the further collapse of the neutron star into a black hole. 
The delay could be of the order of months or years or perhaps as 
short as hours (see Zeh et al. 2004). Nevertheless, the time delay 
here is negative and it is much more plausible that the compar- 
ison of light-curve morphology to SN 1998bw may need more 
than two parameters to be accurate. 

Fits of SN 1998bw templates were also employed to study 
the colour evolution of SN 2010bh in detail as presented in Fig. 
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6, where the curves are shown without the empirical afterglow 
component. From the first two data points, it is clear from the 
g'- r ' and g’~ z'colours that there is a blue component that can- 
not be modelled by the templates and is interpreted as the shock 
breakout in Sect. 3.1. From these two colours, it is also possi- 
ble to see that SN 2010bh becomes red faster than SN 1998bw 
did. The r'- t'and J-H colours remain roughly constant, which 
shows that the changes occur on a broader wavelength scale. The 
standard colour evolution from blue to red is shown by r' - z' 
and r'- J, which at late times evolve bluer and redder than SN 
1998bw templates, respectively. 

4.4. Comparison with SN 2006aj 

We now compare results from optical data by Ferreroetal. 
(2006), who used the same technique and templates to get the 
luminosity and stretch factors for SN 2006aj. They computed 
luminosity factors in the range 0.62-0.76, which are approxi- 
mately in the same range as SN 2010bh without considering the 
measurements in g'z'. Their stretch factors ranged from about 
0.62 to 0.69, which makes the SN 2010bh optical light curves 
wider than those of SN 2006aj. In contrast to the definition of 
the stretch factor, namely that earlier peak times tend to corre- 
spond to a narrower light curve, peak times are earlier in the 
case of SN 2010bh. Whilst the BV photometry for SN 2006aj 
peaked roughly 9 and 1 1 d after the burst, respectively, the g'r' 
photometry for SN 2010bh peaks approximately 8 and 9 d after 
trigger, respectively. This supports the statement that SN 2010bh 
has evolved more rapidly than any other GRB-SNe, given that 
SN 2006aj had the fastest evolution until the discovery of SN 
20l0bh. The JH light curves of SN 2006aj (Cobb et al. 2006; 
Kocevski et al. 2007) are scaled to the luminosity distance of 
SN 2010bh for comparison. After the host-extinction correction, 
SN 2010bh turns out to be as bright as SN 2006aj in the NIR 
as well. Hence, SN 2010bh is similar to SN 2006aj in terms of 
light-curve shape and luminosity. 

5. Bolometric light curve 

The bolometric light curves of SNe are an essential tool for ex- 
amining global luminosity features and enable us to compare 
with other SNe and theoretical models. However, it is difficult 
to obtain such a light curve because of the limited information 
at UV, infrared, and radio wavelengths. Only a gt/aa-bolometric 
1 ight curve can be constructed by using a broad spectral coverage 
to derive a total flux that is then used as a proxy of the bolomet- 
ric flux. To accomplish this task, we employed the wavelength 
range covered by our g'r’i'z'JH filters, i.e.. from 350 to 1800 
nra. 

In the case of SN 2010bh, there are no UV constraints. In the 
following description, no attempts at correcting for the UV flux 
were made. By using the afterglow-subtracted g’r’i'z’JH pho- 
tometry corrected for host-galaxy extinction, monochromatic 
fluxes for each bandpass were derived. Sets of three bandpasses 
were defined to interpolate their corresponding monochromatic 
fluxes using the Simpson's rule. The second-degree polynomial 
result of the interpolation was then integrated over frequency in 
the range of each set of bandpasses. Finally, the total flux in 


the range from 350 to 1800 nm was determined by adding up 
the integrated fluxes of all sets of bandpasses. The total flux is 
transformed to a quasi-bolometric luminosity using a distance 
of 240 ± 17 Mpc to the host galaxy of SN 2010bh (see Sect. 3.1 
for more details). No attempts of extrapolation beyond the lim- 
its of the g' and the H bandpasses were made. Corrections for 
the NIR flux at late times were found to be the most significant. 
The data at to + 30.7 d (rest frame) were corrected for the //-band 
non-detection by assuming that the fraction of H- band flux com- 
pared to the total bolometric flux remains constant at 8% starting 
from to +24. 1 d. Similar corrections were performed for JH non- 
detections at 38.6 and 78.8 d after the burst under the assumption 
of a constant JH flux fraction of 27% at to + 30.7 d. Results are 
shown in Fig. 7 along with quasi-bolometric light curves from 
SE and other GRB-SNe. 

The analysis of the quasi-bolometric light-curve morphology 
yields a peak luminosity of 4.3xl0 42 erg s" 1 at about 8 d after 
the trigger (equivalent to M bo i = -17.87), i.e., approximately 
two times fainter and six days sooner than for SN 1998bw. 
Our luminosity is 16% higher than that computed by Cano et al. 
(201 lb), although consistent to within our 11% of uncertainty. 
Whilst the early peak of SN 2010bh correlates with its narrow- 
ness and low luminosity, this is the case for neither the entire 
GRB-SNe sample nor the local sample of SE SNe (Zeh et al. 
2004; Richardson et al. 2006; Richardson 2009). The morphol- 
ogy of the light curve is similar to that of SN 2006aj (Pian et al. 
2006), although 21% fainter. The peak time also resembles that 
of the type-Ic SN 19941 (Richmond et al. 1996), although SN 
2010bh has a much wider light curve, which is 77% brighter at 
maximum. In terms of peak luminosity, SN 2010bh is similar to 
the broad-lined Ic SN 2009bb (Pignata et al. 2011). It also un- 
derwent the most dramatic late-time decay in the sample, which 
implies that its envelope became rapidly optically thin to y-rays. 
The last statement is supported by the extremely high expansion 
velocities measured for SN 2010bh of the order of 30,000 km 
s~‘ (Chomock et al. 2010a). Another clear feature is the sudden 
decrease in luminosity at around to + 30 d, which contrasts with 
the smooth decay in the comparison SNe at similar stages. This 
indicates either that the atmosphere becomes rapidly optically 
thin or that the assumption of a constant NIR contribution after 
~ to + 31 d underestimates the flux in the JH bands. 

5.1. Physical parameters of the explosion 

We followed the approach described in Valenti ct al. (2008) 
to derive the physical quantities that characterise the explo- 
sion, i.e., we modelled the early and late light curves sepa- 
rately. The early-time phase corresponds to the photospheric 
regime for which the analytical model developed by Arnett 
(1982) has been adopted, initially used for SNe la and adapted 
to SE SNe (e.g., Taubcnberger et al. 2006; Valenti et al. 2008; 
Pignata etal. 2011; Benettietal. 2011). At late stages, the at- 
mosphere becomes nebular, i.e., optically thin, and the emitted 
luminosity is powered by the energy deposition of: (1) -y-rays 
from 36 Co decay, (2) y-rays from electron-positron annihilation, 
and (3) the kinetic energy of the positrons (see appendix A in 
Valenti et al. 2008). However. Maeda et al. (2003) noted that the 
two-component configuration leads to inconsistencies between 
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Fig. 7. Quasi-bolometric light curve of 
SN 2010bh produced by using GROND 
g'r'i'z' JH filters (black open circles) 
in the rest frame. For clarity, the uncer- 
tainties in the host-galaxy extinction 
are not included in the error bars. A 
single black continuous line represents 
the best-fit model (see Sect. 5. 1 ). Early 
and late components of the model 
are smoothly joined at to + 30 d. SE 
and other GRB-SNe quasi-bolometric 
light curves have been plotted as a 
comparison sample: SN 20031w (GRB 
031203: Malesanietal. 2004), SN 
1998bw (GRB 980425; Galamaetal. 
1998b), SN 2006aj (GRB 060218; 
Pian et al. 2006), the broad-lined Ic SN 
2009bb (Pignata et al. 201 1 ), the type-Ic 
SN 19941 (Richmond et al. 1996). the 
type Ibc SN 2008D (XRO 080109; 
Modjaz et al. 2009; Soderberg et al. 
2008), and the type-Ic SN 2002ap 
(Gal-Yam et al. 2002; Foley et al. 2003; 
Yoshii etal. 2003). 


the parameters derived from fitting the early and late light curves 
of SNe Ic, caused mainly by the model limitations in varying the 
y-ray trapping over time. To enable low and high y-ray trapping 
at early and late times, respectively, Maeda et al. (2003) divided 
the ejecta into a high-density inner region and a low-density 
outer region. The emission from the outer region dominates the 
total emission in the optically thick regime at early times, and 
that from the inner region, which has a higher y-ray opacity, 
dominates in the nebular phase at late times. Here, we use the 
same procedure to model the g'r’i'z' JH quasi-bolometric light 
curve of SN 2010bh. 

Given the model explained above, a total of four free pa- 
rameters were used to fit the quasi-bolometric light curve of SN 
2010bh: the total mass of 56 Ni produced in the envelope Mnu the 
total ejecta mass M ej , the fraction of mass in the inner component 
I'm, and the fraction of kinetic energy in the inner component /£-. 
The kinetic-encrgy-to-ejected-mass ratio of the outer region was 
fixed by using its correlation with photospheric velocity at peak 
luminosity (Arnett 1982) 


This expression assumes that the density of the ejecta is homo- 
geneous and that the inner component does not contribute to 
the emitted luminosity in the optically-thick regime. Since the 
photospheric velocity was not available directly from observa- 
tions. the velocity measured by fitting P-Cygni line profiles was 
used as a proxy of y p h. However, the envelope layer where the 
blue-shifted absorption line forms does not necessarily coincide 
with the position of the photosphere, as found when measuring 
different expansion velocities from absorption lines of different 
species. The spread can amount to several hundreds km s~* (see 
Jones et al. 2009 for an example of type-II SNe). For the spec- 


tra of SN 2010bh, Chomocketal. (2010a) obtained velocities 
of about 35,000 and 26,000 km s" 1 from the Si n T6355 feature 
roughly 21 and 6 d after the burst, respectively. Since there is 
no measurement at the time of maximum light, which is about 
8-9 d after the trigger, a range of expansion velocities was used 
in the modelling. Assuming that the photosphere lies at deeper 
layers than those where lines are formed and recedes in mass ex- 
posing deeper and slower layers, in the modelling we employed 
photospheric expansion velocities of 2.5, 2.8, and 3.1 xlO 4 km 
s~‘. Other physical and mathematical quantities such as opac- 
ity and integration constants were chosen to be the same as in 
Cano et al. (201 lb). 

Our fitting procedure consisted of two steps. First of all, we 
modelled the data around maximum luminosity (5 < t - to < 
30 d) using the Arnett’s model, assuming that only the outer 
component contributes to the total luminosity at this stage. We 
then obtained A/Ni.out and M e j.out- Secondly, we modelled the late- 
time data (t - % > 30) using the nebular-phase components of 
Valenti et al. (2008) assuming that both the inner and outer re- 
gions contribute to the total emitted luminosity. Here, we fixed 
A4ii,out and A^j.out to the values obtained in the first step of the 
fitting procedure and only fu and fz were allowed to vary. 

From the above modelling scheme, our best-fit parameters 
were M nloal = (0.135 ± 0.001)A4, M ej>out = 2.37-2.90 M 0 , 
fy = 0.36 ± 0.04, and fg = 0. 1 1 - 0. 15 for the three selected ex- 
pansion velocities at the photosphere. All of these results com- 
bined together provided the total masses and energy of the ex- 
plosion listed in Table 3 (statistical errors only). The 56 Ni mass 
is independent of the chosen expansion velocity at the photo- 
sphere, A/s, = (0.21 ± 0.03 )M 3 , given that it is proportional 
to the luminosity. In contrast, since the ejected mass changes 
significantly, the weighted mean and the RMS of the results of 
the three models in Table 3 were employed to compute a fi- 
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Table 3. Physical parameters of the explosion with a varying expansion velocity. 


^ph 

[km s~'] 

^Ni.out 

[A4] 

A^ej.out 

[A4] 

4f 


M Ni 

[Ms] 

My 

IMA 

£k 

[ 1 0 52 erg] 

25.000 

0. 135 ± 0.001 

1.52 ±0.05 

0.36 ± 0.02 

0.15 ±0.03 

0.21 ± 0.02 

2.37 ±0.10 

1.85 ±0.09 

28,000 

0.135 ±0.001 

1.70 ±0.05 

0.36 ± 0.04 

0.13 ±0.07 

0.21+0.03 

2.64 ±0.14 

2.52 ±0.22 

31.000 

0.1 35 ±0.001 

1.87 ±0.06 

0.36 ± 0.02 

0.11 ±0.02 

0.21 ±0.02 

2.90 ±0.12 

3.37 ±0.13 


nal value of M ej = (2.60 ± 0.23 )M©. The total kinetic energy 
and energy fraction were derived in the same way, implying that 
E k = (2.4 ± 0.7) x 1 0 52 erg and f E = 0. 1 2 ± 0.02. 

Using comparable independent data, Cano et al. (201 lb) fol- 
lowed a similar procedure to derive physical quantities from 
quasi-bolometric data. Whilst our values for Mj are consistent 
with those in Cano et al, (201 lb) to within 1.5cr, they found that 
Mni = (0.10 ± 0.01 )M© for SN 2010bh, which is two times 
lower than our value. The discrepancy affects the determination 
of the kinetic energy as well, which is connected in Eq. 3 to u ph 
and M e j. The causes of these inconsistencies are: (1) the host- 
galaxy extinction employed, which is E{B - T) hoS t = 0.39 + 0.03 
mag in our case and 0.18 ± 0.08 mag in Cano et al. (201 lb); (2) 
the inclusion of an inner component hidden at optically thick 
stages, which increases the 56 Ni mass by a factor of fu\ and (3) 
the choice of different expansion velocities of the photosphere. 
Furthermore, even when including the uncertainties in the host- 
extinction determination, which are <r,w N ~ 0.02 both here and 
for Cano et al. (201 lb), the discrepancy persisted. Nevertheless, 
evidence for the existence of a dense inner layer of0.94±0. 1 5M?, 
is provided and supported by the sub-luminous post-maximum 
phase of SN 2010bh, which is indicative of a high trapping of 
y-rays. 

The large amount of 56 Ni produced is consistent with the 
class of SN associated with GRBs. but also matches the value de- 
rived for the highly energetic type-Ic SN 2004aw (Valenti et al. 
2008), which in contrast had a much broader light curve that 
varied more slowly with time. Most interestingly, SN 2010bh re- 
sembles SN 2006aj in terms of light-curve shape. The 56 Ni mass 
produced in the explosion is practically the same (Mazzali et al. 
2006), whilst the ejected mass is 20% higher for SN 2010bhand 
the total kinetic energy is significantly different (~ 2 x 10 51 erg 
in the case of SN 2006aj). This ensures that SN 2010bh is re- 
markable in terms of expansion velocity, which has been one of 
the greatest ever measured (Chomock et al. 2010a). 

We note that the end of the optically-thick regime and the 
beginning of the nebular phase cannot be accurately defined. 
The rule of thumb is that the photospheric phase ends 30 d after 
explosion at the earliest and that the nebular phase starts 60 d 
after the explosion at the latest. Since Arnett’s model fits our 
data relatively well until day 30, it is defined as the end of the 
optically-thick era. Nevertheless, we are unable to clearly es- 
tablish whether the period between 30 and 60 d after the burst 
corresponds to the nebular phase already, as the ejecta are most 
probably neither completely thick nor sufficiently thin. 

The model reproduces the luminosity at maximum bright- 
ness and thereafter, despite underpredicting the luminosity dur- 
ing the rising phase. This is because the light curve is domi- 
nated by shock-breakout emission at early stages. The same is 


obtained in Sect. 4, when trying to fit templates of SN 1998bw 
to our multicolour light curve, where an extra component for 
the shock breakout is required. The light curves of SN 2010bh 
peak so early that an even more rapid evolution than observed 
is expected after maximum luminosity. This early-wide peak di- 
chotomy could be explained by the 56 Ni distribution in the enve- 
lope. If there were more 56 Ni produced near the surface, the peak 
would be early and wide, i.e., the SN would rise more slowly but 
much earlier. In contrast, if the 56 Ni were concentrated towards 
the centre, it would lead to a much later peak, although shatper 
rise, producing a narrower shape of the maximum (M. Bersten, 
2010, private communication; Nomoto et al. 2010). 


6, Conclusions 

Spanning a time range from 12 hours to 83 days after the trig- 
ger and covering from 190 to 2300 nm in wavelength (see 
Sect, 2), we have presented UV/optical/NIR photometric data 
of XRF 100316D/SN 2010bh. Given the results introduced and 
discussed in Sect. 3, 4, and 5, we have drawn the following con- 
clusions: 

- Broad-band SEDs at early times demonstrate the existence of 
red and blue components identified as synchrotron emission 
from the XRF afterglow and the cooling envelope after shock 
breakout, respectively. 

- A significant amount of dust along the line of sight through 
the host galaxy (A^ost = 1.2 ± 0.1 mag) is consistent with 
the two-component SED model and agrees with the faint de- 
tections at UV wavelengths. 

- By comparing with earlier X-ray results from Starling et al. 
(2011), we have demonstrated that the temperature of the 
blackbody component decreases with time, which is consis- 
tent with a scenario of a cooling expanding atmosphere. 

- By performing an additional analysis of the thermal compo- 
nent and the earlier X-ray measurements we have measured 
expansion velocities that are consistent with SN expansion 
and an initial apparent emission radius of 7 x 10 u cm. This 
radius is slightly larger than the size of WR stars, which are 
the most likely GRB-SNe progenitors. If a WR star were 
the progenitor of XRF 1003 16D/SN 2010bh, then the initial 
radius could indicate that there was a massive dense stellar 
wind surrounding the progenitor. 

- Our multicolour light curves after host-galaxy correction and 
subtraction of the afterglow component, have peak r'i ' lu- 
minosities of about 0.5 -0.7 times that of SN 1998bw and 
consistent with those of SN 2006aj. Similarly the NIR lumi- 
nosity at maximum is as bright as SN 2006aj and 0.6 -0.7 
times that of SN 1998bw. The excess in the g' band indi- 
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cates that SN 2010bh has a hotter photosphere than that of 
SN 1998bw at the time of maximum brightness. 

- We have found that SN 2010bh is the most rapidly evolv- 
ing GRB-SNe to date, reaching maximum luminosity 8-9 
days after the burst in the g'r' bands. At late times, it also 
fades more rapidly than SN 1998bw showing redder colours 
as well. This behaviour is also evident in the bolometric light 
curve, which decays faster than for any SN in the comparison 
sample. 

- The physical parameters of the explosion are derived by 
means of the quasi-bolometric light curve constructed from 
our g'r'i'z'JH photometry. The modelling is performed us- 
ing Arnett’s model (Arnett 1982) for data around peak and 
standard y- ray deposition at later times. A high-density in- 
ner component with roughly 26% of the total mass is re- 
quired to reproduce the flux ratio between maximum lumi- 
nosity and tail. The total mass of 56 Ni produced in the en- 
velope is Mn i = 0.21 ± 0.03 M s , which precisely matches 
the value derived for SN 2006aj, whilst the total ejecta mass 
of Me j = 2.6 ± 0.2 M q exceeds the value for SN 2006aj by 
20%. However, the kinetic energy turns out to be higher at 
E k = (2.4 ± 0.7) x 10 52 erg, making SN 2010bh the second 
most energetic GRB-SN after SN 1998bw. 

The association between XRF 100316D and SN 2010bh is 
particularly interesting, since for the second time the cooling of 
the shock breakout has been detected in a GRB-SN. It is also 
unique in revealing a hot component that possibly contributes 
even at blue maximum brightness, one of the largest host-galaxy 
extinctions measured for this kind of transient, and the fastest 
rise among GRB-connectcd SNe. 
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Appendix A: Optical and near-infrared photometry 


Table A.l. GROND photometry of field stars used for relative photometry. 


R.A.(J2000) 

r h . m . s j 

Dec.(J2000) 

g' 

r' 

i r 

z' 

J 

H 

K, 

07:10:30.10 

-56:14:58.1 

17.981(16) 

17.423(10) 

17.246(12) 

17.065(16) 

16.587(12) 

16.737(13) 


07:10:34.37 

-56:15:58.1 

18.092(16) 

17.365(10) 

17.158(11) 

16.973(16) 

16.477(12) 

16.646(13) 

16.964(31) 

07:10:28.77 

-56:14:30.9 

16.221(13) 

15.831(10) 

15.742(11) 

15.596(15) 

15.224(07) 

15.404(08) 

15.933(15) 

07:10:31.17 

-56:14:45.8 

16.232(12) 

15.535(09) 

15.305(10) 

15.094(15) 

14.606(06) 

14.642(07) 

15.083(11) 

07:10:29.97 

-56:16:24.4 

17.703(13) 

16.435(10) 

15.902(10) 

15.570(15) 

14.875(06) 

14.864(07) 

15.216(11) 

07:10:30.98 

-56:14:26.3 

16.804(13) 

16.391(10) 

16.267(11) 

16.117(16) 

15.733(09) 

15.949(09) 

16.374(22) 


Notes. All magnitudes are in the AB system and corrected for neither Galactic nor host-galaxy extinction. Errors are statistical only. 


Table A.2. GROND photometry of XRF 100316D/SN 2010bh after image subtraction. 


Time Interval g' r' r" z' j H 

[d] after the trigger 


0.49023-0.49470 

0.56046-0.60240 

1.47634-1.49123 

1.69356-1.70396 

2.49804-2.51041 

3.50557-3.52164 

4.48027-4.48963 

5.46324-5.47229 

7.47679-7.48174 

10.47338-10.48259 

11.47550-11.47994 

14.47912-14.49638 

15.48101-15.48997 

19.52885-19.53123 

21.48630-21.49543 

25.49771-25.51874 

32.48587-32.50697 

40.46892-40.50042 

83.43058-83.45161 


20.75(04) 20.90(04) 

20.81(03) 20.91(03) 

20.75(04) 20.90(04) 

20.92(04) 20.77(03) 

20.70(04) 20.58(04) 

20.67(03) 20.40(03) 

20.43(04) 20.21(03) 

20.33(04) 20.11(04) 

20.26(06) 19.96(04) 

20.38(03) 19.91(03) 

20.50(03) 20.01(03) 

20.71(04) 20.07(03) 

20.72(05) 20.07(04) 

21.33(10) 20.36(05) 

21.52(06) 20.48(04) 

21.91(05) 20.78(03) 

22.78(10) 21.42(04) 

> 23.11 21.83(06) 

> 23.95 23.66(23) 


20.87(05) 20.77(07) 

20.94(04) 20.76(04) 

20.59(05) 20.51(05) 

20.69(05) 20.65(06) 

20.50(05) 20.20(04) 

20.24(03) 19.94(03) 

20.15(04) 19.78(04) 

20.04(04) 19.61(04) 

19.87(06) 19.58(06) 

19.98(04) 19.37(03) 

19.99(04) 19.33(03) 

20.08(03) 19.38(03) 

20.16(05) 19.32(04) 

20.39(07) 19.51(05) 

20.48(04) 19.65(03) 

20.66(03) 19.87(03) 

21.28(05) 20.39(04) 

21.74(07) 20.81(05) 

> 22.91 2178(25) 


20.69( 14) > 20.04 

20.71(26) > 20.67 

20.47(22) 20.69(33) 

20.40(32) > 20.09 

19.99(17) >20.09 

19.72(17) 19.94(21) 

20.00(19) 20.22(23) 

19.93(19) 19.92(23) 

19.69(19) 20.20(28) 

19.53(14) 19.85(25) 

19.49(14) 19.94(18) 

19.51(12) 19.67(22) 

19.59(13) 19.85(19) 

19.70(19) > 19.71 

19.64(13) 19.95(19) 

19.67(13) 20.30(21) 

20.07(19) > 20.41 

> 20.34 > 20.43 

> 20.95 > 20.90 


Notes. The template image for subtraction was taken 234 days after burst. The time reference is k = 55271.531 13 MJD, the date of the GRB 
ti'igger. All magnitudes are in the AB system and corrected for neither Galactic nor host-galaxy extinction. Upper limits are all at the 3tr 
confidence level. Errors include the systematics of the photometric calibration. 



